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Abstract 



The paper presents the conservative dynamics of two-body point-mass systems up to the 
third post-Newtonian order (1/c^). The two-body dynamics is given in terms of a higher 
^ ' order ADM Hamilton function which resuhs from a third post-Newtonian Routh functional 

, for the total field-plus-matter system. The applied regularization procedures, together with 

making use of distributional differentiation of homogeneous functions, give unique results for 
the terms in the Hamilton function apart from the coefficient of the term [vpidifr^^ . The 
result suggests an invalidation of the binary point-mass model at the third post-Newtonian 
\ order. 

O : PACS number(s): 04.25.Nx, 04.30.Db, 97.60.Jd, 97.60.Lf 

W)' 1 Introduction and summary 



The calculation of general relativistic equations of motion for compact binary systems, in the 
recent past, was exclusively devoted to the obtention of higher order post-Newtonian gravi- 
ty \ tational radiation reaction contributions. With the works by Iyer and Will [Q], Blanchet 

and the authors the radiation reaction levels have been completed up to 3.5 post-Newtonian 
(3.5PN) order, i.e. to the order beyond the Newtonian dynamics (for the 2.5PN order 

see, e.g., the review [^). Quite recently, Gopakumar et al. |5| succeeded in giving first results for 
the gravitational radiation reaction in compact binary systems at 4.5PN order, applying balance 
equations between far-zone fluxes and near-zone losses of energy and angular momentum. The 
conservative 3PN and 4PN orders of approximation, however, have not been tackled so far (up 
to the 2PN level of approximation see, e.g., the review and the paper Q). 

It is the aim of the present paper to develop the two-body point-mass dynamics up to the 
3PN order of approximation using the canonical formalism of Arnowitt-Deser-Misner (ADM) . 
There are several aspects which make such a calculation rather interesting: (i) the well-known 
problem of applicability of Dirac delta distributions for the source of the general relativistic 
gravitational field; (ii) the need of the 3PN dynamics for a better undertanding of the innermost 
stable circular orbit for binary systems |9|, |ll|; (iii) the importance of the 3PN dynamics 
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to control the viability of the 2PN filters for gavitational wave measurements from inspiralling 
compact binaries |12]. 



Within om' canonical approach, together with applying Hadamard's "partie finie" technic, 
some specific analytic regularization formula and a generalization of it achieved by us — the 
non-generalized formula has been successfully used by Damour [^], Damour and Schafer 
and Kopeikin at the 2PN level and by us at the 3.5PN level — , and the distributional 
differentiation of homogeneous functions, we were able to explicitly calculate and uniquely regu- 
larize all terms which occur at the 3PN order of approximation, but one (notice: this is the only 
ambiguous term up to the 3.5PN order of approximation, inclusively). The coefficient of the 
following term turns out to be finite but ambiguous: —{vpidi)^r~^, or, in non-reduced- variable 
form (apart from a factor 4): 

miv\v{ o ^ r, ( Gm2\ , m2v\v'^ ^ ^ ( Gm 
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where mi,v\,rsi and m2,V2,rs2 denote the masses, velocities, and Schwarzschild radii {vsa — 
Gma/2(?, in isotropic coordinates) of the bodies 1 and 2, respectively, and ri2 their relative 
distance. This term describes the quadrupole-quadrupole interaction of the kinetic energy tensor 
of each body with the gravitational tidal field of the other body scaled to the Schwarzschild radius 
of the former one. The bodies, in this interaction term, are obviously not tidally deformed (this 
is expected from Ref. |^ where it is shown that the tidal deformation comes in at 5PN only) but 
they seem to have been attributed extensions of Schwarzschild radius size. Those extensions are 
beyond the applicability of the 3PN approximation so the obtained result may suggest that a 
fully consistent 3PN compact binary model needs extended bodies as source for the gravitational 
field. Only within the complete theory the bodies' Schwarzschild radii are expected to enter the 
equations for the binary system in a consistent and unique manner. 

The paper is organized as follows. In Section 2 we introduce the point-mass model and 
the post-Newtonian approximation scheme and we develop the constraint equations to the 3PN 
order of approximation. To obtain an autonomous Hamilton function for the bodies at 3PN 
order, dropping the dissipative 2.5PN level, we introduce in Section 3 the Routh functional 
for the total field-plus-matter system and eliminate the field degrees of freedom for the bodies' 
variables. The autonomous Hamilton function comes out of higher order. The Section 4 is 
devoted to the explicit calculation of the higher order autonomous Hamilton function applying 
the regularization procedures of the Appendices B.l and B.2. In the Section 5 a thourough 
investigation of the obtained ambiguity is undertaken. In the Section 6 we compare our results 
with limiting expressions known from the literature. The Appendix A presents several explicit 
metric coefficients and some useful formulae for inverse Laplacians. The Appendix B.l is devoted 
to the Hadamard's "partie finie" regularization. The Appendix B.2 presents a powerful analytic 
regularization procedure based on an analytic formula derived by Riesz and generalized by us 
in this paper. In the Appendix B.3 the analytic regularization procedure devised by Riesz is 
given. The Appendix B.4 shows the distributional differentiation of homogeneous functions. 

We use units in which IGvrG = c = 1, where G is the Newtonian gravitational constant and c 
the velocity of light. We employ the following notation: x = (x*) {i = 1, 2, 3) denotes a point in 
the 3-dimensional Euclidean space 'S? endowed with a standard Euclidean metric and a scalar 
product (denoted by a dot). Letters a and b are body labels, so G denotes the position 
of the ath point mass. We also define Fq := x — x^, := jraj, := ra/va-, and for a ^ b, 
Tab ■= Xq — Xf), Tab '■= |rafe|) ^ab '■= ^ab/^ab] \ ' \ stauds here for the length of a vector. The linear 
momentum vector of the ath body is denoted by pa = (pai), and denotes its mass parameter. 
Indices with round brackets, like in 4>(2)j give the order of the object in inverse powers of the 
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velocity of light, in this case, We abbreviate (5 (x — Xq) by 6a- An overdot, like in x^, 

means the total time derivative. The partial differentiation with respect to x* is denoted by di 
or by a comma, i.e., di(j) = cp^f, the partial differentiation with respect to we denote by dai- 



Throughout this paper we extensively used the computer algebra system Mathematica [16|. 



2 The constraint equations up to 3PN order 

We consider a many-body point-mass system which interacts with the gravitational field ac- 
cording to the theory of general relativity. For such a system the constraint equations in the 
canonical formalism of ADM read (see, e.g., Eqs. (2.8) in [17U) 
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Here gij = ^g^j are the field variables (the prefix denotes a four-dimensional quantity, all 
unmarked quantities are understood as three-dimensional), g^^ is the inverse of gij {g^-' gjk = ^l), 
g := det{gij}, "|" indicates the covariant derivative with respect to g^j, R is the curvature scalar 
formed from the metric gij] vr*-' is the canonical conjugate to the field gij. Spatial indices are 
raised and lowered using g^^ and gij, respectively. 

We use the following coordinate conditions (see, e.g., Eqs. (2.9) in |]j 
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where hj^^ is the transverse traceless part of gij — 6ij. The trace- free field momentum vr*-' can 
be split into two parts, a longitudinal tt*-' and a transversal tt*-'"'""'" one: 
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where vr*-' can be expressed in terms of a single vector vr* as follows: 



TT 



= vr^ + vrf. - 6ij7r% + A ^11%,,. 
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If both the constraint equations (|l|)-(§) and the coordinate conditions (P)-(@) are satisfied, 
the Hamiltonian of the theory can be put into its reduced form, which can be written as 



H 
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The reduced Hamiltonian contains the full information for the dynamical evolution of the canon- 



ical field and matter variables M, 18 



We expand the constraint equations (|^) and (^) in powers of 1/c where we take into account 
that (see, e.g., Ref. Q) 
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To calculate the reduced Hamiltonian (0) up to 3PN order we have to expand the Hamiltonian 
constraint equation (||) up to Making use of Eqs. @, (Q), and (|8|), after long calculations, 

we obtain 
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We also need to expand the momentum constraint equations (|2|) up to l/c' . Using Eqs. (^), 
), and (P), we get 
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All functions entering the right-hand sides of Eqs. (|9|) and ( [lO| ) can be written as sums of 
terms of different orders in 1/c. To the orders needed in our calculations, they read 
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Using Eqs. (ll)-([l^) we extract from Eq. (^) the Hamiltonian constraint equations valid at 
individual even orders in 1/c, up to the order 1/c^. They are 

-A0(2) = ^ma6a, (15) 
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"2^ V^(4)iij + 4 i^(4)»i,fej + 2^(4)ii,fc^(4)ifc,i - 2^(2),rjh(^6)ij- 



(18) 



Explicit solutions of the equations (|T5|) , (|T6|) , and ( p^ ) for the functions (/)(2) , <A(4) > aiid (/)(g) , 
respectively, are shown in Appendix A. The full solution of the equation (|l8|) for the function 
0(8) is not known. We split ^(g) into two parts 

4>{8) = (PiS)! + 4>{8)2- (19) 

The function 0(8)i is explicitly calculable and can be found in Appendix A. The unknown part 
of (/)(8) we call (p(^s)2- It is given by 

4'{8)2 = ^~"^5'(8)2) 

where 
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Using the Eqs. (|ll|)-(|l^ we extract from Eq. (p^ the momentum constraint equations valid 
up to the order 1/c^: 
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where 
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2)^j5) + (^'^(2) + '^(4)) ^Js) + <^(2)^;5r + 2vrJ)^5];fc 

(25) 



Explicit solutions of the equations (|2^ ) and ( p3[ ) for the functions tt^^^ and 7r|g^ , respectively, 
are given in Appendix A, where one can also find the formula for the function vr^g^ connected 

with vf by means of Eq. (^) . 

The TT-part of the metric in the leading order fulfils the equation (see, e.g., Eq. (14) of [^) 

= ^I™A(^4)ki, (26) 
where ^ 

a °' 

and where Sjlj^^^ is the TT-projection operator defined by (see, e.g., Eq. (4) of Q) 

-{6ki - A-'dkdi){6ij - A-'didj) ] . (28) 



TTkl 



The explicit formula for the function hj^-j can be found in Appendix A. 

In our calculations we also need the explicit formula for that part of the metric function 
^(6)ij which diverges linearly at infinity. So we write 

_ > TTdiv , > TTconv /r,qx 

where hj^^g^ ~ r and /iJI™''" ~ 

as r ^ oo. The function hj^^" is given by the integral 

(see, e.g., Eq. (12) in [|) 

. TTdiv / .N 1 srTTfc/ f^3^r ^^^(4)fc; / , a i / i rTTH f A^'^ A \ +\ /"?n^ 

h(fi)ij (x, t) = Jdx — — — [x,t) \x-x \- d^j ('^)^' ) ( ' ' ( 

where the function ^(4)ij is defined in Eq. (27). We have calulated the integral from the right- 
hand side of Eq. (^). The result is given in Appendix A. 



3 The conservative Hamilton function up to 3PN order 

In the following we are interested in a conservative 3PN Hamiltonian which depends on body 
variables only. To achieve this goal we firstly transform our field-plus-matter Hamilton func- 
tional (0) into a Routh functional which is a Hamilton function for the bodies but a Lagrange 
functional for the field (notice for the following the crucial difference between a Hamiltonian 
and a Lagrangian: the functional derivative of the latter is zero, the one of the former not). The 
Routh functional up to 3PN order turns out to be 



-R<3PN 



Xa , Pa , hj^,-j , hj^ij ] = i7<3PN - y" d\ vrg^^ /i JJ.^. , (31) 
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where i/<3PN is the Hamiltonian up to 3PN order. The equations of motion of the point masses 
determined by the Routh functional (^Tj) read 

5i?<3PN . 9i?<3PN 



Pa 



dyia 



dpa 



(32) 



We ehminate now in the Routh functional (^) hj^-j and hj^-^ by Xq, Xq, p^, and Pa through 
solving the field equations which result from ff<3PN (see Eq. (p6|)). After that the Routh 
functional ( |3l| ) becomes a higher order matter Hamilton function (denoted by a tilde) of the 
variables x^, x^, Pa, and pa'. 
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Xaj Pa? h(4:)ij (^a) Pa) 



TT 

(4)ii 



(Xa? Pa; Xqi, Pa) 



The equations of motion of the bodies determined by the Hamilton function (33) read 

0^<3PN _ ^ ( dH<sP^ \ 

dpa dt \ dpa ) ' 



Pa 



^-g<3PN ^ d ( gg<3PN \ ^ 
9Xa dt I 9Xa y ' 



(33) 



(34) 



In these equations the higher time derivatives may be eliminated by applying lower order equa- 
tions of motion. The elimination of the higher derivatives in the Hamilton function ( |3^ ) would 
result in a redefinition of the body variables (see, e.g., Ref. pl|]). 

The reduced Hamiltonian i^opN of Eq. ( ^T| ) can be obtained from Eq. (^). After dropping 
full divergences (including the term ^'-'vr*-'"^'^, which can be written as [(2-7r-' + A~^7r^.^)7r*-'"'"^]^j, 
cf. Eq. (^)), it reads 
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Le us note that the integrals of the full divergences we have dropped in Eq. ( |35| ) do not contribute 
to the Hamiltonian because they fall off at infinity at least as l/r"^. This can be inferred from 
the following asymptotic behaviour of the functions entering on the right-hand side of Eq. (P) 
(cf. Ref. p 
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see also our discussion below Eq. 

The higher order Hamiltonian ff<3PN of Eq. (^) can be split as follows 

at Pat Pa ) Pa ) 

+if2PN (Xa, Pa) + -f^SPN (x a J Pa J Xfl) Paj ; 



(37) 



where the Hamiltonians Hq through -ffspN can be extracted from Eqs. ( ^T[ ) and ( |35| ) by means 
of the expansions (|l3)-(|l|). The Hamiltonians Hq through H2PN are known. We re-calculate 
them below for completeness. Our aim, however, is to calculate the 3PN Hamilton function 

-f^SPN- 

The calculation of the Hamiltonians Hq through //2PN can be performed directly as we 
explicitly know all functions needed to perform the integrations. The Hamiltonians can be 
written in the form 



Ho = Jd\^ma5a, 
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The Hamiltonian -ffspN, as extracted from Eqs. ( [31] ) and (^), reads 
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(42) 
(43) 
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(45) 
(46) 
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1 J. I.TT tTT PaiPaj , TT e,/^ \i' 
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For the calculation of the Hamiltonian H-^p^ we perform some manipulations which allow us 
to do the integrations in Eq. (^) without explicit knowledge of all functions entering on the 
right-hand side of Eq. (p7[). There are three such functions: the part 0(8)2 of the function i;^(8) 
[cf. Eq. (p^j, the function vf^^^, and the part hj^^"^"^ of the function hj^^j. 
The function 0(8)2 we eliminate by means of the identity 
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To eliminate the unknown function nlL we use the following relation, which can be proved 
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All the terms which depend on the unknown function hj^-j we were be able to write as a 
full divergence [see Eq. ( |60| ) below]. To do this we have used the following relation, which can 
be derived using the explicit formula (^) for the operator and the traceless property of 

the function hj^^.y. 

{^™A(4)ki - ^(4)ii) hj^ij 

= { ^ (A"^^(4)A:fc,j) - 2 (A~M(4)jfcj,^ + i (^A"2A(4)fc;jfc;) ^J^ij} j ' (^^) 
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In the last stage we eliminate the field momentum vr^i^"'" by Xq, x^, Pa, and pa- We use the 
field equation for the field momentum in the leading order. It reads (see, e.g., Eq. (13) of 1^]): 

One can check that the following relation holds 

('^(2)^(3)) = '^(2)^(i) + 2^5) • (52) 
Substituting Eq. (||) into Eq. (||) we obtain 



vr 



Sr = I + I '^(2) ^(3) + ^(5) • (53) 

Using Eqs. (p^), (pl|), (^), (^9[), (^), and ( [5^ ) we rewrite the Hamilton function -ff3PN given 
by Eq. (|47| ) as a sum of terms 

6 _ 
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7=1 
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4 Results of regularization procedures 

To diminish the number of terms we perform calculations in the center-of-mass reference frame, 
so we can use the relations 

Pi + P2 = 0, pi+p2 = 0. (61) 

We also use, at the 3PN level, the Newtonian relations between the coordinate velocities xi, X2 
of the bodies and their relative velocity v := xi — X2: 



xi 



m-2^ 



X2 



mi 
'Id 
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where M := mi + 771-2 is the total mass of the system. 

To shorten the formulae we introduce the following reduced variables 
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r Pi P2 
n := -, p := — = , q 



Pi 
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where 
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We also introduce the reduced Hamilton function 



(63) 
(64) 

(65) 



H depends on masses of the binary system only through the parameter v. From now on the hat 
will indicate division by the reduced mass ^. 

To calculate the Hamiltonians Hq through -?/2PN we proceed as follows. The terms Hq, H^, 
Hii, and H21 [given by Eqs. (|38D, (|39|), (p|), and (44), respectively] we regularize by means 
of the Hadamard's procedure described in Appendix B.l. The terms H12 and 7^22 [given by 
Eqs. ( ]4^ ) and (|45|), respectively] we regularize using the procedure from Appendix B.2. The 
term 7^23 [from Eq. (46)] is a full divergence. The integrand in H23 falls off at infinity as 
1/r^ so H23 does not contribute to the Hamiltonian. We have checked that the regularization 
procedure of Appendix B.2 applied to H23 gives zero. The final result for the Hamiltonian 
Hf) + i^N + -^iPN + -ff2PN coincides with that known in the literature (see, e.g., Eq. (3.1) in [^). 

The calculation of the Hamiltonian -ffspN is much more complicated. The term H31 of Eq. 
(|55|) we regularize using the Hadamard's procedure described in Appendix B.l. After long 
calculations we obtain 

^31 = (-5 + 351/ - + 35z.3^ (p2)4 

+1| (-7 + 42Z. - 53u' - 6u') {p'f + (1 - [2(n • pf{p'f + 3(n • pW] }i 
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^1 (-108 + 515z/ + 84z^2) (p2)2 + (46 + 35i/)z^(n • pfp'^ - (45 - 255z/)z/(n • P)^|^ 

47:1 f-500 + 17601/ - 2317/^^) - 3(480 - 1897z/)z^(n • p)44 + 4:(2 - (66) 
.60 t V / ) r-^ 16 



The calculation of ^^32, ^^33, and H^^ [defined by Eqs. (pq)-(]58|)] we do by means of the 
regularization procedure described in Appendix B.2. The performance of the procedure differs 
considerably in the case of H^2 ^-nd H^^ compared to H^^. In and -f^ss each term can be 
regularized separately. This is not the case for i/34. For each term of taken separately the 
limit (130) of Appendix B.2 does not exist. By series expansion with respect to e one can see 
that for each term entering H34 the sum from Eq. ( |130| ) contains a part proportional to 1/e 
(similar divergences arise in the dimensional regularization procedure in quantum field theory, 
see, e.g., [^). After collecting all 1/e terms (they are of the type ([75|) only) the parts cancel 
each other. Furthermore, after regularization, there are terms of the type (^) proportional to 
lnri2. After collecting all such terms the logarithms cancel out, too. Unfortunately, the result 
of application of the regularization procedure to is not unique. The ambiguity which arises 
can be expressed in terms of exactly one unknown number and will be discussed in detail in the 
next section. 

The result of applying the regularization procedure of Appendix B.2 to Eqs. (|56[)-(p^), after 
tedious calculations, is 



H: 



32 



3 

64 



z/ |(-4 + 113i/)(p2)2 + (12 - 31iy){n ■ p) 



H V 

512 



6037r2 - 9632 - 4 (l004 + 157:^) 1/ 



P 



+3 



4640 - 6037r2 + 12 (^5tt^ - 44) u\ (n • p)H ^ + 



1 1 



(67) 



H: 



33 



^ |4(n • p)(n • q)(p • q) - 5(n • 



5(n • p)2q2 - (n • p)^(n • q)^ + 13p^q^ + 2(p • q 



2„2 



1 2 



n • v)(p • q) 5p^ + (n • p) 



(npj 



p2(q-v)+2(p-q)(p-v) 



-3(11 • p)^(q • v) + (n • q) [p^ + (n • pf\ [(n • p)(n • v) - (p • v) 

+ 4r f (4(17 - 10u){n • p)(n • q)(n • v) - 2(15 + 22u){n ■ pf{n ■ q) 
48 t 

-2(51 - 8z/)(n • q)p2 - 4(6 - 5z/)(n • p)(p • q) - 4(1 - 2i/)(n • v)(p • q) 



4(7 - 2i^) (n • p)(q • v) + (n • q)(p • v) 



8(n • p)(n • v)p2(p • v) 



+8(n • pf{n ■ v)(p • v) + 2(n • pfp^^ + 5{pYv^ + (n • p)V - 5(n • ^r)\pY 
-6(n • p)2(n • v)2p2 - 5(n • p)\n ■ v)^ - 4(n • pf{q_ ■ v)^ - 4p2(q • v)^^ ^ ^ 

1^ |5(5 - 7i/)(n • pfin ■ v) + 10(3 - 5i/)(n • pfin ■ vf 
+3(17 - 35i/)(n • p)(n • v)p2 - 28(3 - 8z/)(n • p)(n • v)(p • v) 
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-15(2 - 3z^)(n • p)2(p . v) + 2(24 - 77i/)(n • pfv^ + 2(9 - 29z^)(n • ^^fp^ 
-3(4 - 9z^)p2(p • v) - 2(12 - 37i/)p2v2 + 4(6 - 17z^)(p • ^f^^ 



+ 



1 



1536 
6 



ui3 



mi? - 10832 36(48 - hTx^)v (n • p)(n • v) 

2\ o/'/ic:„2 AaA\,, („ ,,^2 , ^^ac\^\ nT7-^2 



+2 



3(16 vr^) - 2(457r^ - 464)i/ (n • v)^ ^ 11600 - 9277r^ 20(97r^ - 80)z^ 
176 37r^ 6(176 - 157r^)i/ " ' ^ 



34 



+ 



|9(19 - 13i/)(p2)2 - 18(3 - 43z.)(n • p) V - (5 - 303z.)(n • p)4|l 
-180(103 evr^) (33944 2257r2)z^ 



1 



1920 



P 



-F27 



60(27r2 - 3) - (2216 257r2)z^ (n • p)^ 
p2-3(n.p)2 



96 



(926 - 637r2)i/ 



1 



(69) 



The number <jJ in Eq. ( |69D is unknown due to the non-uniqueness of the results of the regular- 
ization procedure, as discussed in the next section. 

The integrals given by ^^35 and i^ae [defined in Eqs. ( |59| ) and (|^)] are integrals of full 
divergences. To discuss them let us first observe that the 3PN Hamiltonian given by Eqs. (|54|)- 
( pOD not only describes the system of two point masses but can be applied also to the system 
of two dusty bodies (provided the following substitutions are made: m-a^a Pa, Pai^ 
(p^/ma)5a ^a/Pa, etc, where pa and Pa are the mass and the linear momentum density of 
the ath body, respectively). It is obvious that for dusty bodies the integrands in H^^ and H^q 
are locally integrable, so these terms can contribute to the Hamiltonian only if the integrands 
fall off at infinity slower than 1/r^. 

Now we proceed as follows. We check whether the integrands in and H^q fall off at 
infinity at least as 1/r^. If yes, the terms do not contribute to the Hamiltonian, if no, we must 
study them in more detail. Because the asymptotic behaviour of all functions entering and 
H3Q is the same for both point masses and dusty bodies, we use the above prescription also for 
the system of point masses, treating them as a limiting case of extended dusty bodies. 

One can check that the whole integrand in ( |59|) falls off at infinity as 1/r^, so formally 
ifss does not contribute to the Hamiltonian. For checking the consistency of our regularization 
procedures we have integrated out all terms in H^^ except the last two ones (which contain the 

and (/>(8)2)- To do this we have used the regularization of Appendix B.2. 

This nonzero result is connected with 
the ambiguity of the term as explained in the next section. 

The integrand in H^q does not fall off at infinity fast enough not to possibly contribute to the 



unknown functions vr^^-j 
We have obtained zeros for all terms but the first one. 



Hamiltonian. The reason is that the function hTJ:. . is a sum hj^'^" + hTJ:^°^'^ , 



where h 



TTdiv 



~ r 



1/r as r ^ 00 (cf. Eqs. ( ]29| ) and (^0[)). The integrand in Hi 



and 

falls off at infinity as 1/r^ but it falls off only as 1/r xui 



36 



?,TTdi\ 



for hll^ 



^TTconv 



Therefore we have 



calculated H^q for h 



TT 

"(6)^ 



I, TTdiv 



' for h^., 

We have used the regularization procedures from Appendixes 



(4)U 



B.l and B.2 (the need to use the Hadamard's procedure of Appendix B.l arises because Ah^-"^ 
contains some terms with Dirac deltas). To calculate AhJ^-j properly we also had to employ 
the rule of differentiation of homogeneous functions from the Appendix B.4. The result is zero. 
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On the basis of the above discussion we put the integrals given by H^^ and H^q equal to 
zero, 

^35 = ^36 = 0, (70) 

and adjust the regularized expression different from zero to the Hamiltonian 

Collecting the Eqs. (|6^)-(|70|) we finally obtain the autonomous higher order 3PN Hamilton 
function ifspN- It reads 



^3PN(r,p,v,q) = -L (_5 + _ 70Z.2 + 35j,3^ (p2) 



+ 42u - 53iy'' - Qu-" ) {p'y + (1 - 2i^)iy'' 2(n • p)^(p^)^ + 3(n • p)*p 



1 



48 



<{ 3 ( -27 + 140z^ + 96i/^ ) (p^)^ + 6(8 + 25i/)z^(n • p)^p^ - (35 - 267iy)u{n • p)M ^ 



2„2 



+ 



1 



1536 



-4800 - 3(8944 - 3157r2)z^ - 71361^^] p^ + 9 (2672 - 31571^ + 224i/) u{n • P)^ 



1 



+ — <{ 12 + (884 - 637r')iy }>- + i^33 (r, P, v, q) + a; p^ - 3(n • p)^ 



(71) 



where the number u! is unknown and the function H^^ (r, p, v, q) is given by the right-hand side 
of Eq. (H). 



5 Ambiguity 

The source of ambiguity given by the regularization procedure described in Appendix B.2 can be 
explained as follows. Via integration by parts, some terms in H34 can be represented in different 
ways. The regularization method applied to both representations give different results. 

As an example let us consider the integral of 4'{2)4'{2),i4'{2),jhJ^-j which, by means of integra- 
tion by parts, can be written as integral of — ^</'(2)'^(2),ii^^')ij (™ Appendix B.4 another such 
term is treated). Application of the regularization procedure to the difference between these 
integrals (the integrand is a full divergence) gives (in the reduced variables) 



I Jd'x (^0?2)0(2),,/iS,) ^ = f [P' - 3(n • P)'] ^- (72) 



Only if the result ( |7^ would be zero, application of the regularization procedure to the integrals 
of 4>{2)4>{2),i4>{2),jh'^. and -^<A(2)'A(2),ji^(4)ij would give the same results. 

We can also obtain the result (^) in a different way. Let us denote the integrand from the 
left-hand side of ([72|) by Fi^i and let us consider the volume integral 



'd^Fi^i, (73) 

B(0,i?)\[B(xi,ei)UB(x2,e2)] 

where B{-Xa,ea) (a = 1,2) is a ball of radius Sa around the position Xa of the ath body and 
B (0, R) is a ball of radius R centered at the origin of the coordinate system. We apply Gauss's 
theorem to the integral (^) and then we calculate the limit 

lim (p dai Fi + lim <f> dat Fi + lim i) dai Fi (74) 

ei-*0 / £2^0 / R^oo J 

aB(xi,£i) aB(x2,£2) dB{0,R) 
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with the normal vectors pointing inwards the spheres dB (xa,£a) and outwards the sphere 
dB{0,R). It is easy to check that the integral over the sphere dB(0,R) vanishes in the limit 
i? — > oo whereas the integrals over the spheres dB {xa, £a) diverge as Sa — 0, so to calculate 
them we use the Hadamard's procedure from Appendix B.l. The result coincides with that 
given on the right-hand side of Eq. ([7^). 

We have checked that for all these terms entering for which one can use different, via 
integration by parts, representations the ambiguity is always a multiple of the quantity 



,2 



p^-3(n-p)^ =-{upAy-, (75) 



r 



or, written in non-reduced form and with the momenta substituted through the velocities (apart 
from a factor 4), 

where r^i and rs2 denote the Schwarschild radii (r^a = Gmajlc? ^ in isotropic coordinates) of 
the bodies 1 and 2, respectively. This expression indicates that the radius of the bodies might 



come into play already at 3PN as indicated e.g. in Ref. |2C]. The interaction described by the 
terms (^) is the interaction of the (Newtonian) kinetic energy tensor of each body with the 
(Newtonian) tidal potential of the other body scaled to the respective Schwarzschild radius. 

We have tried to resolve the ambiguity (^) in several ways. Firstly we have extracted the 
local non-integrabilities in i?34. Let us denote by F the total integrand of the B.^/^ and by F^a 
(a = 1, 2) the non-integrable part (i.e., of the order 1/r^ or higher) of the Laurent expansion of 
the F around the position Xq of the ath point mass. Then we replace F by the locally integrable 
expression F — {Fs\ +-^s2)- We have checked that the results of applying the regularization 
procedure of Appendix B.2 are the same for F and F — {Fs\ + -^52) (so the regularized value 
of the singular part F^x -|- of the integrand is zero). This is so because by subtracting the 
local non-integrabilities we have transferred the integrand, initially locally non-integrable but 
integrable at infinity, to an integrand locally integrable but non-integrable at infinity. 

We have studied two further possibilities to overcome the ambiguity. The first one, based 
on the Riesz's kernel representation of the Dirac delta distribution, is described in detail in 
Appendix B.3. The Riesz's kernel regularization of the divergence ( |72|) gives zero, i.e. it takes 
into account also the contributions coming from the points x^. We were yet not able to compute 
all terms in 7^34 using the Riesz's kernel regularization because of serious calculational problems, 
as described in Appendix B.3 (only the first and the momenta dependent part of the fifth term 
out of six terms in ff34 we could calculate, but with ill-defined (1/e and lnri2) multiplication 
factors for the terms of type (|75|)). In the second case we have tried to employ the rule of 
differentiation of homogeneous functions coming from the distribution theory. This is described 
in Appendix B.4. Here we also haven't fully succeeded in removing the ambiguity. 

We conjecture that the ambiguity has its origin in the zero extension of the bodies. We 
started with point-like bodies but the formalism reacted in such a manner that the Schwarzschild 
radii of the bodies got introduced. They, however, are far beyond the applicability of the post- 
Newtonian approximation scheme and thus, the result turned out to be ambiguous. In the future 
we shall investigate our conjecture in more detail. 

We close this section with another observation. According to the applied canonical formalism 
one is allowed to put K^'^ and tt'-'"'""'" equal to zero at any freely chosen initial time. The energy 
of the initial state is given by /f<3PN(xa, p^, /i^"'^ = 0,7r*-''^^ = 0). It exactly contains the first 
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term of But this term is ill-defined even with respect to the Riesz's kernel regularization 
procedure (Appendix B.3). This is a further indication that the binary point-mass model is 
ill-defined at the 3PN level. 



6 Comparison with the known results 

The 3PN Hamiltonian derived in Section 4 we now compare with the results known in the 
literature. We have found only two such results: the test body limit and the static part of the 
full Hamiltonian. Let us stress that these comparisons do not fix the ambiguity present in our 
paper. 

The Hamiltonian (expressed in the reduced variables defined in the beginning of Section 4) 
describing a test body orbiting around a Schwarzschild black hole reads [0] (we restored the 
explicit dependence on the speed of light c) 



2c^r / V 2c^r / V V 2c^r / 



(77) 



Expanding ( [77[ ) with respect to 1/c and taking the 1/c^ contribution we obtain the 3PN Hamil- 
tonian for a test body: 

fvtest _ Hp')' 7ipy 27{pr 25p2 1 

^3PN- ^ + ^- 

To obtain the test-body limit of the 3PN Hamiltonian of our two-body point-mass system we 
have to substitute u = into the right-hand side of Eq. ( |7l| ) (note that the ambiguous term and 
the term ^^33 vanish for u = 0). The result coincides with (|78|) . 

The static part of the Hamiltonian is defined by means of the conditions 

Pa = 0, Tr'^^^ = 0, hj^ = 0. (79) 

In the paper |23] one can find the following expression for the static part of the many-body 
point-mass 3PN Hamiltonian: 

a h c d e 

X j ■■■ j d^Xi . . . d^X5 6 (Xi - Xa) S (X2 - Xfe) 6 (X3 - Xc) 6 (X4 - Xd) S (X5 - Xe) 
/ 2 



+ + , (80) 

Vn2?'23?'34?^45 n2f'23f'34?^35 n2?'23?'24r25 / 

where r^n := — x„|. We have calculated H^^^^ from (^) for a two-body system. For 
this we have used the Hadamard's regularization procedure of Appendix B.l. The result is (in 
reduced variables) 

To calculate the static part of the 3PN Hamiltonian in our approach, we have to implement the 
conditions (79). After that the quantities H^2^H2,& from Eqs. (|5^)~^(|60[) vanish. Only some part 
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of the expression H31 of Eq. ( |55|) survives. Applying for this part the regularization procedure 
from Appendix B.l one obtains the formula (|8ll). 
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A Some explicit solutions 

In this appendix we cite the known approximate solutions of the constraint equations (|l]) and 



(|2|). The functions 0(2), 0(4), 0(6) (implicitly), Trjgw and hJJ-- can be found e.g. in Refs. [^] and 



|24|. The functions 0(8)i, vf^^j, and hj^^" have been calculated by us. 

In most of the Poisson equations given below there are source terms of the form J2a /(x)5a, 
where the function / is singular at x = Xq. The Poisson integrals for these terms we have 
calculated as follows 



|5:/(X)5,| = Aeg(xJ5,| = ^/,eg(xJA"l5, = -i-^/,,g(xjl 



(82) 



where the regularized value freg of the function / is defined in Eq. ( Ill ) of Appendix B.l. 
The Poisson integrals which do not contain Dirac deltas we have computed using the explicit 
formulae (|103| )~ (|108| ) for inverse Laplacians listed at the end of this appendix. 



The Hamiltonian constraint equation (|T|) is explicitly solved up to the 1/c^ order. The 



solutions of Eqs. (15) and ([iq ) for the functions 0(2) and 0(4) are known for n-body point-mass 
system. They read 

a " 

^ - 2 rngmb , 1 V- Pg 



The solution of Eq. (^7j) for the function 0(g) is fully known only for two-body point-mass 
systems. It can be written as 

4>{6) = '^'(6)1 + 4>i6)2 + (85) 



where 

4>{6)1 



(86) 
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Computation of the Poisson integrals from Eqs. (|86| ) and (^) gives the results 



+ 



327r ^ m^ra (167r) „ ^ "^"^a'' 



(89) 



V ) a b^a 



+ [-8 (Pa • Pfe) (Va • Vb) + 3 (pa • (pfo • Vfo) - 8 (p^ • Vfo) (pfc • V„)] A 



.-1 



ran 



+4 (p, • V,) (p5 • V,) (V, • V5) A-^^ 

n 



(90) 



The Poisson integral needed to calculate (/>(6)3 for two-body point-mass systems reads (here 



A 



-1 



1 mainb / 4 



{3rt - Urlrab + 18r>2^ - 12r„r3, + 3r 



4 



32 (167r)2 rSr^r^^ 
+28rj^rb - 12r^rafcrfc - Uvarl^n + 28r^{,rb - 30r^r^ OOrar^fcr^ 
-30r^br^ - 3Qrarf - SQrabrf + 35r^ 
3 1 



647r mar^ 



p; - 3 (ria - Pa) + 



1 1 

327r nib 



2(p.-V,)^^-4^ 



+ 



2pf (V • Vb)2 - 4 (pfe • V) {pb ■ Vb) (V • Vb) 



-4 (Pfe • V) (pb ■ Vb) (V • V5)] ( A-i^ ) + 8 (p, • V)^ ( A 



2 / A-l 



ran 



+ g(Pb-Vb)2(V-Vb)M A- 



-I'b 



Only part of the solution of Eq. (O) for the function is known, namely 



where 



5, 



(8)1 



V- I 1 .3 1 . . 1 . ^ 5 ^ 15 . \ p2 



(91) 



(92) 
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-^9(2) A — — + —ri/ 



64"^"^ m4 16 m« '2 ml ^^'^'^ j 
We have computed the Poisson integrals from Eq. (p^). They are given by 



(8)1 



1 



E 



2\3 



+ 



1 



EE 



(P^)' r 1 



+ — 



(93) 



1 ^ (Pa)^ + ^ • P")'' Pa + 3 ("afe " P") 



rriambrabra 



+ E E ^ (pa f- - -) + 15 (n„, • pj^ f- + -H 



+ 



(16vr) „ ^_^^ v„ft 



1 



EE 



1a ^{0(4),,/iS,} 



1 >p "^a"T-b 



The momentum constraint equations (|22D-(P4|) can be written in the form 



(94) 



(95) 



Making use of the decomposition (^) it is not difficult to write the solution of Eq. ( p5| ) in the 
form 

1 



TT ■ 



A; ' 



where the operator is defined as follows 



-^5ijdk + dikdj + Sjkdi - ^didjdkA ^ ) A ^ 



(96) 



(97) 



Using Eq. (^) we have obtained the solutions of Eqs. (|2^) and (^) valid for n-body point- 
mass systems. The solution of Eq. (22) reads 



4 = i^E^^<^^{-'^^.(;^ 



+ 2 



(98) 



The function vr^g^ , connected with the function vrjg^ by means of Eq. (^) , is equal to 



1 



(3) 
1 1 



TT 



(3) - Svr 



(99) 
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The explicit formula for the function vr^g^ reads 



TT, 



(5) 



~T [P'^i^a + PajK) - ("a ' Pa) ^ij + {^^a ' Pa) 



(16^)^ V ^a 



+ 



(16vr) a b^a [ \raJ,kn 



(Pa • 



+4 (p, • V) {didaj + djdai) + 4 {pajd, + Pa^^,) (V • Vj] A 



.-1 



+ 



^Si, {pa ■ Va) (V • V„)' - {didaj + ^,^a^) (Pa • V,) (V • V,) 



n 



-Adid, (pa • V) (V • Va) [ J-") + Id^dj (p, • (V • v„)' f A-2^ 



. (100) 



The TT-part of the metric in the leading order (1/c^) is the solution of Eq. (|2i 
by (here Sab '■= + fb + Vab', be aware that in Eq. (20) in ^ there is a misprint) 



It is given 



h 



TT 

{4)ii 



1 1 1 r 

iI^E;;^{[Pa-5(na-Pa)^]5.,+2p, 



+ 



3{na ■ Paf - 5pl nln^ + 6{na ■ Pa)(n>aj + n^Pai)] 



o (167r) a b^a ^ ^ 



+2 



Ta^Tb , 12 



' ah 



+ 



1 /r.^ 



I ab' a I nb \ a 



1 1 



+ 



Vn / r, 



^ab 



17 ^ 4 _^8/l_^4N^ 



'^ab^a '^a^b ^ab \^a '^abJ 



(101) 



The part of the TT- metric at 1/c^ which diverges linearly at infinity is given by Eq. (pO|). 
After computing the double Poisson integral in (|30|) one finds 



LTTdiv 



24 167r 



^ — I 7(na • Pa)^ - llp^ 5ij + 26paiPaj 



7Pa - ("a • Pa)^ - 10(na • Pa)(n>ai + n^l^ai)} 



1 fl 



1 / 3 2 \ 
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+— {didbj + djdbi - didaj - djdai) — + , , , o ^i^j 
12 rnf, 144r?, 



Below we list explicit formulae for the inverse Laplacians which have been used throughout 



this appendix: 




A- 1 = 

Tan 


In Sah, 


A- 1 = 

Tan 


36 V 


n 




n 


7200 ^ 



3 



(103) 



12 



+90rab4 + 63r,4) + ^ (^a + 2r>2, _ 3r^, + Irlr^ + Gr^.r^ - 3r,4) Ins.b, (106) 
A-V.rfe = ^ (63r^ + OOr^r^^ - 62r2r2^ - 30r,r3, - r^, + 90rln - QOrlrabH + 30r,r2,rb 
-SOrlbU - 126ryb - dOraVabrf - Q2rlbrl + OOr^rj^ + OOr^fer^^ ^ 53^4 
+ t4t (-3^a + 2rlrl, + r^, + Gr^r^ + 2r2,r2 - 3r,4) In s^b. (107) 



120 



= (-63r^ - 90r;Jr,, - 2r>^, - 90r,r;J, - 63r:, + 150r>, 

+90r>aferfe + DOrar^^rb + Ih^^rlbU + I26rlrl + dOrarabvl 
+I26rlbrl - 90r,rf - gOr^^rf - 63rfe^ 

+ ^ (3r^ + 2rlrlb + 3r^, - Gr^r^ _ gr^.r^ + 3r,4) In Sab- (108) 

B Regularization procedures 

In this appendix we describe techniques which can be used to regularize integrals which appear 
in our paper. More details are given in Ref. ||25||. 
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B.l Hadamard's "partie finie" regularization 

Let / be a real valued function defined in a neighbourhood of a point Xq G K^, excluding this 
point. At Xo the function / is assumed to be singular. We define the family of auxiliary complex 
valued functions /n (labelled by unit vectors n) in the following way: 



/n : C 9 e ^ /n(e) := / (x^ + en) G C. 
We expand /„ into a Laurent series around e = 0: 

oo 
m=-N 



(109) 



(110) 



The coefficients am in this expansion depend on the unit vector n. We define the regularized 
value of the function / at Xo as the coefficient of in the expansion ( |11C| ) averaged over all 
directions: ^ 

/reg(xo) := — <l>dnao{n). (Ill) 

We use formula ( p.llD to compute all integrals which contain Dirac delta distribution. It means 
that we define 

y'A/(x)5(x-x,) := Aeg(x,). (112) 

The procedure described here was used by one of the authors (cf. Appendix B in [|l^) in 
the calculation of the 2PN and 2.5PN ADM Hamiltonians for many-body point-mass systems. 
More details on the applications of the Hadamard's regularization can be found in [25|. Also 
the relation of the Hadamard's procedure to the regularization described in Appendix B.2 and 
the rule of differentiation of homogeneous functions of Appendix B.4 is discussed in [p5| . 

B.2 Riesz's formula based regularization 

The following formula, ffistly derived by Riesz (see Eqs. (7) and (10) in Chapter 2 of ||2^), can 
serve as a tool to regularize a class of divergent integrals of rf using the analytic continuation 
arguments: 

r(^] r(^) r (-^±^\ 
' \ ^ ) \ ^ ) \ 2 ; a+P+Z 



d x rf r: 



3/2; 



rcg 



12 



(113) 



In deriving the 3PN Hamiltonian for a two-body point-mass system more general integrals 
than those which can be regularized by means of the formula ( |113| ) appear. We succeeded in 
deriving the generalization of the formula (|113| ) for those integrals. The generalized formula 
reads 



J d\ rf (ri + r2 + ri2) 



R{a,P,-f)r 



a+/3+7+3 
12 ' 



(114) 



where 



R{a,l3,j) 



2-K 



-a 



/? - 7 - 4) 



^(-7) 

/1/2 (a + 2, -a - 7 - 2) + h/2 (/? + 2, -/? 



7 



-I1/2 (q + /3 + 4, -a - /? - 7 - 4) - 1 



(115) 
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The function in Eq. ( [LISD is defined as follows 



A/2 {x,y) ■■ 



_ Bi/2 {x,y) 



B{x,y) 



(116) 



where B stands for the beta function (Euler's integral of the first kind) and -B1/2 is the incomplete 
beta function; it can be expressed in terms of the Gauss hypergeometric function 2F1: 



1 



Bi/2 (x, y) = 77— 2Fi[ 1 -y,x;x + 1;- 



2^x 



1 



(117) 



The regularization procedure based on the formulae ( |114| ) and ( |115D consists of several steps. 
We enumerate them now. 

1. By means of Eqs. (^) and ( |62D we replace P2 by — pi, P2 by — Pi, and xi and X2 by v. 

2. We eliminate the unit vector n2 by the relation 

(118) 

'2 

Then we expand the scalar product ni • 1112: 

(119) 



ri , ri2 
n2 = — ni H ni2. 

r2 r2 



ni • ni2 



2riri2 



3. All integrands we reduce to integrands which depend on ri and r2 only. To do this we use 
a set of substitutions. For the integrands of the type (ni • pi)'^ f {f 11^2)1 k = 1, . . . , 6, we 
use the rules (obtained by considering the integration in prolate spheroidal coordinates) 



(ni •pi)/(ri,r2) a/(ri,r2), 



where 



(ni • pi) /(ri,r2) 
(ni • Pif f{ri,r2) 
(ni • Pi)^/(ri,r2) 
(ni • Pif f{ri,r2) 
(ni •pi)V(n,r2) 

a := (ni2-pi) 



a^ + h^]f{n,r2) 



a^ + ^ab^'^f{n,r2), 

a' + 5aH' + h'y{n,r2), 

a5 + 5aV + ^a5^)/(ri,r2), 



a'+^-^a%' + ^aH' + ^b^]fin,r2), 



ro — r 



1 ~ ' 12 



2riri2 



(120) 
(121) 

(122) 

(123) 

(124) 

(125) 

(126) 
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After this step the integrand becomes the Unear combination of the type 

N 



Y.cir'^'r^'{n+r2+n2r' , (128) 



7=1 

where aj, /?/, and 7/ are integers; the constant coefficients cj may depend on ri2, pf, pf, 
v^, (ni2 • pi), (ni2 • pi), (ni2 • v), (v • pi), (v • pi), and (pi • pi). 

4. We define the auxihary function 

Jria, /3, 7) := R{a + ^e, /? + ,^e, 7) r-^Y'3+.+{^+.)e+3 ^^29) 
and for the combination (|12§| ) we calculate the limit 

limY,ciJriai,l3ini)- (130) 

It turns out that for all terms considered in the paper this limit is of the form 

A + B^ + C-, (131) 

where A, B, C do not depend on fj, and i^. 

5. As the regularized value of the integral of the sum ( p. 28] ) we take the number A, i.e. 

N 



d^x ^ cir'^'r2' (n + r2 + ri2)'" 
•' 7=1 



:=A. (132) 

reg 



B.3 Riesz's kernel based regularization 

Three-dimensional Dirac delta distribution can be represented as the limit 

5(x-x,) = limr(x,x,), (133) 

where J"(x, x') is the Riesz's kernel defined as follows: 

r(x,x') := ^ %^J x-xr-'- (134) 

2 



The regularization of a class of integrals based on the kernel ( 134 ) and its relation to Hadamard's 
procedure of Appendix B.l is discussed by Riesz |26]. For us, the Riesz's kernel regularization 
consists in using the kernel ( [134D instead of the Dirac delta in the source terms of the constraint 
equations, solving these equations, and performing the regularization of integrals needed to 
obtain the 3PN dynamics in the way described below. 

For a point-mass two-body system we need two families of Riesz's kernels, each associated 
with the different mass. Let the first family be labelled by the parameter a, and the second one 
by f3. After replacing the Dirac deltas by Riesz's kernels the right-hand sides of the constraint 
equations (|l])~(|2|) depend on a and /3. The first step of the procedure is to solve constraint 
equations iteratively keeping a and P as unspecified complex parameters. 
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This is a very difficult task. Out of all equations for the functions (/>(„) we have fully solved 
only the simplest one for the function <p{2)- Also only the leading order momentum constraint 
equation for vf^g.^ we have solved. The TT-part of the metric is not known even in the lowest 
order — only the part depending quadratically on the momenta has been calculated. 



The solutions obtained by us are enough to calculate that part of the full divergence (72) 
(discussed in Section 5) which depends quadratically on the momenta. The regularization of the 
full divergence by means of the procedure described in Appendix B.2 gives a non-zero result. If 
one uses the Riesz's kernel generalization of the function (j)i^2) ^"^^ of the momenta dependent 
part of the function hj^^ - the result of regularization of the full divergence using the procedure 
described below is zero. 

The procedure is similar to that described in Appendix B.2. The first three steps of the reg- 
ularization procedure from Appendix B.2 apply here. After that all integrands to be regularized 
are functions of ri, r2 and the parameters a and (5. So any integrand can be written as 

/(ri,r2,a,/3). (135) 



We need now to use formulae analogous to the generalized Riesz's formula ( |113| ) from Appendix 
B.2 to perform integration of (135) for any a and /3. Let's denote the result of integration by 



R{a,P). Then the following limit is calculated: 

lim R{a + /xe, /? -|- ue) . (136) 

For the divergence (|72|), to perform the integration it is enough to use the formula (|113|) from 



Appendix B.2 for 7 = 0. The limit ( 136 ) comes out to be zero. 



For all terms for which we were able to apply the above procedure we obtained the limits 



(|136| ) of the form given in Eq. (131). However, we were not able to check the cancellation of all 
divergent terms in i?34 after regularization because we were not able to calculate all terms using 
the Riesz's kernel procedure. But again, the ill-defined terms we have obtained were of the type 



B.4 The rule to differentiate homogeneous functions. 

In distribution theory there exists a rule to differentiate some homogeneous and locally non- 
integrable functions under the integral sign [^^. We have studied how useful this rule is in our 
calculations. 

Let / be a real valued function defined in a neighbourhood of the origin in M'^. / is said to 
be a positively homogeneous function of degree A, if for any number a > 

/(ax) = aV(x). (137) 

Let k := —A — 2. If A is an integer and if A < —2 (i.e. A; is a non-negative integer), then the 
partial derivative of / with respect to the coordinate has to be calculated by means of the 
formula (cf. Eq. (5.15) in ^) 

where {dif)^^.^ means the derivative computed using the standard rules of differentiations, S is 
any smooth close surface surrounding the origin and dai is the surface element on S. 
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As an example of applying the rule (|l3q ) let us consider the full divergence (connected with 
the fourth term in the i/34 part of the 3PN Hamiltonian, see Eq. (|58D) 



Jd\ I (2vrj3) + A-ivr[3) (2^3) + A-ivr[^)_, . /.JT^j ^ . (139) 
Applying the regularization procedure from Appendix B.2 to the divergence (|139D gives a result 



much more complicated than that given by Eq. (75). After performing in ( |13S| ) differentiation 
with respect to one finds that only in one term the rule (|13^) can be applied. The term reads 



The rule (138) applied to Trf^-. yields 



3 (27ri3) + A"V[3)^,,) /iS,vrf3)^ (140) 
(3),ifc 

Jd\ {-i (2vr|3) + A-V[3),,,) (141) 

The integral ( |l4lD is calculated by means of Hadamard's procedure from Appendix B.l. After 
adding the result of regularization of the integral (|141| ) to the result of regularization of the 
integral ( |139| ) we obtain a multiple of the quantity (175|). 

Using similar considerations we were always able to restrict the ambiguity to a multiple of 
the quantity (|75|). Let us also stress that the part of the 3PN Hamiltonian is written in 
such a form that there is no need to use the rule ( |138| ) for individual derivatives appearing in 
H34. 

A method of applying the rule ( |141| ) in regularizations of integrals can be found in Section 



6 of |25|. We have yet found that this method (without modifications) is not able to give zero 



for all full divergences one meets in our calculations. 
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